On a C. de Boor's Conjecture in a Particular Case and 

Related Perturbation 



Zhe Li, Shugong Zhang*, Tian Dong 

School of Mathemathics, Key Lab. of Symbolic Computation and Knowledge Engineering 
{Ministry of Education), Jilin University, Changchun 130012, PR China 



Abstract 

In this paper, we focus on two classes of D-invariant polynomial sub- 
spaces. The first is a classical type, while the second is a new class. With ma- 
trix computation, we prove that every ideal projector with each D-invariant 
subspace belonging to either the first class or the second is the pointwise 
limit of Lagrange projectors. This verifies a particular case of a C. de Boor's 
conjecture asserting that every complex ideal projector is the pointwise limit 
of Lagrange projectors. Specifically, we provide the concrete perturbation 
procedure for ideal projectors of this type. 
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1. Introduction 

Polynomial interpolation is to construct a polynomial p belonging to a 
finite-dimensional subspace of ¥[x\ from a set of data that agrees with a given 
function / at the data set, where F[a:;] := F[xi, . . . , x^] denotes the polynomial 
ring in d variables over the field F. It's important to make the comment that F 
is the real field M or the complex field C in this paper. Univariate polynomial 
interpolation has a well developed theory, while the multivariate one is very 
problematic since a multivariate interpolation polynomial is determined not 
only by the cardinal but also by the geometry of the data set, cf. [H, 0]. 
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Recently, more and more people are getting interested in ideal interpola- 
tion, which is defined by an ideal projector on F[a;], namely a linear idem- 
potent operator on ¥[x] whose kernel is an ideal, cf. [3|. When the kernel 
of an ideal projector P is the vanishing ideal of certain finite set S in F*^, 
P is a Lagrange projector which provides the Lagrange interpolation on H. 
Obviously, P is finite-rank since its range is a dimensional subspace of 
¥[x]. Lagrange projectors are the standard examples of ideal projectors. 

It's well-known that an ideal projector can be characterized completely 
by the range of its dual projector, cf. 0, [^, @, 01- Given a finite-rank linear 
projector P on ¥[x], the kernel of P is an ideal if and only if the range of its 
dual projector is of the form 

with some finite point set H C F'^, D-invariant finite- dimensional polynomial 
subspace C ¥[x] for each ^ G H. 6^ denotes the evaluation functional at 
the point ^, and Q^[D) will be explained in the next section. 

In the univariate case, for an integer n, there is only one D-invariant poly- 
nomial subspace of degree less than n, which implies that every univariate 
ideal projector can be viewed as a limiting case of Lagrange projector, cf. 
[sl. This prompted C. de Boor to define Hermite projector as the pointwise 
limit of Lagrange projectors and pose the following conjecture in [9|. Indeed, 



this conjecture had been raised in [10|] with certain restriction. 



C. de Boor's conjecture A finite-rank linear projector on C[a;] is an ideal 
projector if and only if it is the pointwise limit of Lagrange projectors. 



B. Shekhtman constructed a counterexample to this conjecture for ev- 
ery d>2). In the same paper, B. Shekthman also showed that the conjecture 
is true for bivariate complex projectors with the help of Fogarty Theorem 
(see [13] )• Later, using the fact that any pair of commuting matrices can 
be approximated by pairs of diagonalizable commuting complex matrices (see 



]), C. de Boor and B. Shekhtman [15| reproved the same result. In 
addition, by Theorem 8 of 16|], 15|] also proved that certain low-rank mul- 



tivariate ideal projectors are the limit of Lagrange projectors. Specifically, 



B. Shekhtman [17[ completely analyzed the bivariate ideal projectors which 
are onto the space of polynomials of degree less than n over real or complex 
field, and verified the conjecture in this particular case. 
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Since for every d > 3, there exist ideal projectors that are not the point- 
wise hmits of Lagrange projectors, the question that what kind of ideal pro- 
jectors can be perturbed as Lagrange projectors lies ahead. For this purpose, 



B. Shekhtman [18| theoretically presented a symbolic algorithm which can 



determine whether an ideal projector is the limit of Lagrange projectors or 
not. However, as mentioned by this paper, such a method isn't yet feasible 
in practice. 

In the converse case, one wonders how to generate a sequence of Lagrange 
projectors practically such that this sequence of Lagrange projectors con- 
verges pointwise to a given Hermite projector. More deeply, B. Shekthman 



raised the following question in [19 



Problem Let P be an Hermite projector, and Ph a sequence of Lagrange 
projectors such that Ph converges pointwise to P as h tends to zero. Then, 
what's the relationship between the trajectories of the points in varieties of 
kerPh and PI 



In this paper, we deal with two classes of D-invariant subspaces. The 
first one is classical, which has been investigated by many literatures such 



as 



10, 20, 21, 22, 23 



^ The second is some special type, which is inspired by 

examples in 10, J 19 1. We construct a group of interpolation point sets cor- 
responding to each class, and establish the relationship between evaluation 
functionals induced by the point sets and derivative functionals related to the 
corresponding D-invariant subspaces. Based on these results, we generate a 
sequence of Lagrange projectors which converges pointwise to the ideal pro- 
jector with related D-invariant subspaces belonging to either the first or the 
second class. Equivalently, we constructively verify C. de Boor's conjecture 
for the ideal projector of this type. 

The remainder of this paper is organized as follows. The next section is 
devoted as a preparation for the paper. Section 3 and 4 discuss the first class 
of Z)-invariant subspaces and the second, respectively. Finally, main theorem 
of this paper is given in Section 5. 



2. Preliminaries 

In this section, we will firstly introduce some notations and review some 
basic facts about ideal projector. For more details, we refer the reader to 
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Throughout the paper, we use N to stand for the monoid of nonnegative 
integers and boldface type for tuples with their entries denoted by the same 
letter with subscripts, for example, ex. = [ai, . . . , ad)- For arbitrary a G N^, 
we define ||c»:||i = ai + . . . + ad and ck! = ai! . . . a^!- 

For arbitrary tuples a, /3 G N'^, we write f3 < a if ex. — f3 has only 
nonnegative entries. In other words, < is the usual product order on N'^. A 
subset D C N'^ is called a lower set (alternatively, down set, order ideal, etc.) 
if it is closed under <, that is, o; G D implies /3 G D for all f3 < ex. 

Let P be a finite-rank ideal projector on F[a:;]. The range and kernel of 
P will be described as 

ranP ■=P{¥[x]) = {g e ¥[x] : g = Pf for some / G ¥[x]}, 
kerP :=ran(/ - P) = {g e ¥[x] : Pg = 0}, 

where ranP forms a finite-dimensional subspace and kerP a zero- dimensional 
ideal of ¥[x]. It is obvious that the ideal kerP complements the subspace 
ranP , i.e., kerP © ranP = ¥[x]. 

Furthermore, as an infinite dimensional F- vector space, ¥[x] has an alge- 
braic dual F[a3]'. An ideal projector P on ¥[x] also has a dual projector P' 
on ¥[x]', and 



In fact, ranP' is the set of interpolation conditions matched by P. It is easy 
to see that the maximum number of linearly independent functionals in ranP' 
equals the dimension of ranP, namely dimranP = dimranP'. In addition, 
a q = (gi, . . . ,qs) is an F-basis for ranP and A = (Ai, . . . , A^) an F-basis for 
ranP', then their Gram matrix 



is invertible. 

For q; G N'^, we write x^ for the monomial x^^ . . . x'^'^, especially define 
x"" = 1 for a = (0, . . . , 0). Thus, a polynomial p in F[a;] can be expressed as 



ranP' = (kerP) 



± 



{A G ¥[x]' : kerP C kerA}. 



X^q := (Xiqj) 



l<i,j<s 




(1) 
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with p{cx.) G F nonzero. For a polynomial p{x) given as in ([T]), we write 
p{D) : F[a;] ¥[x] 



dx^' ■ ■ ■ dxf 



for the associated differential operator. For a finite-dimensional polynomial 
subspace Q C F[a;], we define 

Q{D) := spanp{g(Z}) : g G Q}. 

We call a polynomial subspace Q a D -invariant polynomial subspace if it's 
closed under differentiation, i.e., for every q G Q, g G Q for alll < z < d. 

Finally, let us recall some facts about combinatorics. For r, m G N with 
< r < m, is the binomial coefficient, i.e., = .. . For a = 

— — 'Vr/ ''Vr/ r\(m—r)\ 

(«!, . . . , ad) G N"^, (^l'"")^,^) signifies multinomial coefficient, that is. 



ai,a2, . . . oli 

Lemma 1. p. 90] Letn,m be arbitrary nonnegative integers satisfying 
n < m. Then 



m—r 



r=0 



f^X n \ ml, n 

r = 



0, < n < m. 



3. The first class of Z)-invariant subspaces 

In this section, we are concerned with the classical D-invariant subspaces 
which only involve directional derivatives. Based on the study of 10|, Isl, |22 



we describe that how the interpolation points are arranged in the straight 
directions, this sequence of evaluation functionals can converge to the corre- 
sponding directional derivative functional. 

Let a. = («!, . . . , ad) G N"^ and p = (pi, p2, • • • , Pd) G (F'^)'^ in which 

Pi = (Pl,l> • • • > Pl,d)' P2 = (P2,l> • • • 5 P2,d), ...,Pd= (Pd,l, • • • , Pd,d) e ^'^ 

are F-linearly independent unit vectors. Then we define a differential opera- 
tor Dp : ¥[x] ¥[x] by the formula 

f)\H\i f 

D'^f = i 

dpT---dpT 
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Furthermore, the next lemma tells that for what type of polynomial, its 
associated differential operator is exactly D^. 

Lemma 2. Let ol = (ai, . . . , ad) G N'^, p = (pi, p2, • • • , Pd) ^ (F'^)'' he as 
above, and let 

d 
i=l 

where ■ denotes the Euclidean product. Then p{D) = D^. 

Proof. The proof can be easily completed by direct computation. □ 

Proposition 3. Given a lower set C N'^ and p = (pi, p2, ■ ■ ■ , Pd) G (F^)*^ 

as above, then the subspace 

d 

Q = spanj,{JJ(pi ■ x)°" : a G 0} C ¥[x]. 

i=l 

is D -invariant. 

Proof. This proposition is the immediate consequence of the differential 
calculus. □ 

The following proposition makes a connection between difference quotient 
and directional derivative of multivariate polynomial, which is a variant of 



Proposition 9.2 of [10| and Proposition 7.2 of [9[. We also provide a simple 



proof, for completeness. 

Proposition 4. Let p = {pi, p2, . . . , pd) G (F'')'^ be as above, h a non-zero 
number in ¥, and let a G N"^, £ G F'^. Then for arbitrary p G F[a^], 



0</3<ci i=l 

where (^) = Y[t=i (^0' ^""^ remainder 0{h) is a polynomial in h. 
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Proof. For convenience, let 



if = 

0</3< 



<»<a fc=l 



Then we can obtain that for m = 0, . . . , ||q:||i, 

^"V(0) V 7M ^71,. 7., ^X^) 

dh"' /'l:! ■■■Pd4 o llTilli Q^ll7dl|i 



7i||i+...+||7d||i=m 



7i,i • • -71,^, • • -,7^,1) • • • )7d,^ 

with 7i = (7i,i, . . . , 7i,d) e 1 < z < 

In the following, we will simplify the right-hand side of the above equality. 
There are two cases which must be examined. 

Case 1: < m < \\cx.\\i — 1. 

In this case, there must exist some I < j < d such that Yl'i=i1ij < '^i- 
By Lemma 1, it follows for such j, 

. -/3,Y"A«J:ti7. 



Thus, for all < m < ||a||i - 1, = q. 

Case 2: m = ||q:||i. 

In this case, if there exists some 1 < j < d such that Yl'i=i > then 
there must exist another 1 < j' < (i, such that X]i=i 7'i,j' < "^i'- Hence, when 
m = \\cx\\i, 

d aj 



j=l /3j=0 ^^■''^ 

if and only if for all 1 < j < d, Yl'i=i li-j — ^j- Combining this with the fact 



/3,.=0 



■3 
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we can conclude that 

'^"""V(O) II II , v;^ 71,1 71, d 7d,i 

dhiHW = ll^lli' ••• 2^ Pi,i ■■■Pd,i ■■■Pi,d ■■■ 

7i,l+---+7d,i="i 7i,d+---+7ti,£i=ad 



7i,i' • • • ' ld,i J \ 7i,rf, • • • , ld4 J c^xi ... 9a; 
That is, 

In sum, we have deduced that 

rfXO) r 0, < m < - 1; 



dh"^ \ \\a\\i\{D'^'p){$,), m=||a||i, 
which imphes that 

Equahty ([2]) follows directly from the above equality. □ 



4. The second class of £>-invariant subspaces 

In this section, we present a new class of /^-invariant subspaces, which are 
spanned by polynomials with special structure. We also provide a group of 
interpolation point sets corresponding to this class of D-invariant subspaces, 
and discuss that how some interpolation points coalesce in the non-straight 
directions, such sequence of evaluation functionals can converge to the dif- 
ferential functional about this class of D-invariant subspaces. 

Before introducing the new class of D-invariant subspaces, we need to 
settle some notations used throughout this section. 

Let a = (oo, ai, . . . , a„) with n > 1 be an n -|- 1-tuple of positive integers 
satisfying 

ao = 1 and ai > . . . > a„ > 2, (3) 

and let 

Cl = (ci,o, Ci,i, . . . , Ci,„), . . .,Cd= {Cdfl, Cd,l, . . . , Cd,n) € F''^^ 
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satisfying that Ci^q, C2,o, • • • , Q,o aren't all zero. Then it's clear that 

n d 

r{ji, . . . , 7d) = J] aj ^ 7ij (4) 

j=0 i=l 

with 7i = (7i,o, 7i,i, • • • , 7i,n) G defines a map r : (N"+^)'^ ^ N. 

Proposition 5. Let a = (oq, ai, . . . , a„), Cj = (cj^oi Cj^i, . . . , Ci^„), 1 < i < d 
and the map r be as above. Let qn,m, ?ti = 0, 1, . . . , ai be polynomials defined 
by 

„ _ ^1 ---^d J\li\\i J\ld\\i 

qn,m- 2^ 7,!...7,!'^^ •••'^'^ ' 

with 7i = (7i,o, 7i,i, • • • , 7i,n) e N"+S c/' = YTj^Q c]y , and let Q be a polyno- 
mial subspace defined by 

Q = span]ir{g„,m : < m < ai}. 

Then the following hold: 

(i) dimQ = ai + 1, 

(ii) Q is a D -invariant polynomial subspace. 

Proof. To prove (i), we assume that there exist ko, ki, . . . , ka^-i, ka^ G F 
such that 

ai — lQn,ai — l ~l~ ^aiQ'n,ai 0. 

Since there exists some I < io < d such that Ci^^o ^ 0, then x™' must belong to 
the support of g„,m. Using this together with the definition of r, we conclude 
that the degree of is m. 

Specifically, belongs to the support of gn.ai, while for all < m < 
Oi — 1, x^^ can't belong to the support of Therefore /c^i =0, which 

implies that 

koQnfl + kiQn^l + . . . + /Cai-lQ'n.ai-l = 0. 

Arguing for m = ai — 1, ai — 2, . . . , 1, as for m = ai, we get fcai-i = 
0, A;ai-2 = 0, . . . , fci = 0, fco = 0, successively. That is to say, the family of 
polynomials ^nm, < m < ai is F-linearly independent. So dimQ = ai + 1. 
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To prove (ii), we firstly notice that if 1 < m < a„, tlien r(7i, . . . , 7^) = m 
implies that 7jj- = with 1 < i < d,l < j < n. Thus, for all 1 < m < a„, 

71.0 7d,0 

^ _ ^1.0 ■■■^d,0 71.0 7d,o 

7i,o! • • • 7(i,o! 
= ^(ci,oa;i + C2,qX2 + . . . + Cd^oa^d)"". 

Consequently, 

= Ci,0^n,m-1, 1 < m < a„. (5j 

Next, let m be an arbitrary integer satisfying an < m < ai, and s the 
minimum integer between 1 and n such that m — > 0. We claim 



Ci,oq: 



n,m- 



1 + ^ Qj-gn.m-aj, a„ < m < Oi- (6) 



This claim together with equality ([5]) immediately means that Q is D-invariant. 

To prove our claim, we will use induction on the number n. When n = 1, 
our claim can be easily verified. Now, assume that our claim is true for n — 1. 
To prove that it holds for n, let k be the maximum nonnegative integer such 
that httn < m, that is, 

kttn < m,{k + l)a„ > m. (7) 

and k' the maximum nonnegative integer such that k'an < m — 1. From the 
definition of qn,m, we obtain that 

k 71,0 7l,n-l 7d,0 7d,n-l 

E/ "-^1,0 • • ■ Cl,n-1 • ■ • Crf^O • • • '^d,n-l 

\ / ^ i i i T 

, ' . , , 7i,o! • • -Ti,™-!! • • -Td.o! • • •Td.n-i! 

(=0 r((7i,o,.--,7l.n-l),.--.(7d,Ov,7d,n-i))=m-«a„ 

71, n 7d,n 

7l,0 + ---+7l,n-l 7d,o + ---+7d,n-i\ ( 1," ■ ■ ■ d,n ^7i,n ^7d,, 

Xi ... X „ 



That is. 



, Ci n! . . . Cd n^- 

7l,n+---+7d,n=' 



gn,,m = ^ yj-(ci,„a;i + . . . + Crf,„Xd)'g„_i,m-Za„, (8) 



/=0 
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which plays an important role in what follows. At this point, we have two 
cases to consider. 

Case 1: a„ < m < a„_i. 

In this case, s = n and < m — /a„ < On-i for all < / < /c. Hence, 

qn-l,m-la^ = J ^^T^^^^'^^^ + • • • + CdflXd)""-^"" , < I < k. 

By iD, we get 

^ 1 

qn,m = ^ 777 rr(Cl,„Xi + . . . + Cd,nXdy {CifiXi + . . . + CdfiXd)'^ 

So we obtain that 
dxi 

1 

=Q,0 ^ 777 _ , .. {Cl,nXl + . . . + Cd,nXd)\ci^0Xi + . . . + C^.O^d)"* ^ + 

L.yTfl i tCtnJ- 

''I 

Cj,n 77 7T77 ^ 7T(Cl,na;i + . . . + Cd^nXd)^ ^{ClfiXi + . . . + CdflXd)"^ 

(/ - l)!(m - /a„)! 

Cj,0'?n,m— 1 ~l~ Ci^fiC[fi,m—a„' 

Case 2: a^-i < m < ai. 
In this case, < s < n — 1. 

By ([7]) and ([3]), it follows that m — kan — a-n-i < 0. Thus, there exists 
some < Iq < k — 1 such that 

m - loan - a.„_i > 0, m - (/q + l)a„ - a„_i < 0. (9) 

As a result, we have that < m — /a„ < a„_i for all /q + 1 < / < A;, which 
implies that 

Therefore, 

dqn-l,m.-la„ _ J Ci,oQ'n-l,m-l-ia„ , ^0 + 1 < ^ < ^^Q) 



ax,: 10, A;' + 1 < / < A;. 
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For all < / < /q, we have a„_i < m — Ian < cii- Then our inductive 
hypothesis implies that 



dqn- 



l,m—lan 



'n-l,m-l-lan + Cj,jrQ'n-l,m-«a„-aj 5 < / < /qj (H) 



dXi 

J=si 



where Si is the minimal integer between 1 and n—1 such that m—lan—cisi > 0. 
It should be noticed that Sq = s. 

By (fTOj) and (ITTi) . we can deduce that 

dQn,m 



dx,, 

k 



CifiC^^ J^iCi^nXl + . . . + Q^n3^d)'q'n-l,m-l-«a„) + 



1=0 



k ^ 



1=1 ^ 

^(Cl,na;2 + . . . + Cd,nXdf ^ Cijg„_i,^ 

Z=0 ■ j=Si 

It remains to show the last row of the above equality and the right-side 
hand of ([6]) are equal. More precisely, for each sq < j < n — 1, let kj denote 
the maximum nonnegative integer satisfying kjttn < m — aj, that is, 

kjan < m — ttj, {kj + l)a„ > m — aj. (12) 

From (|9]), we know kn-i = lo- 
Due to dH]), we observe that 

^ 1 

qn,m-aj = ^ J^ {c2,nX2 + • • • + Cd,nX d)^ qn-l,m-aj-la„, Sq < j < U - 1. 
1=0 

Furthermore, 

qn,m-aj = ^ J^{Cl,nX2 + • • • + Cd^nXdY qn-l,m-aj-la„, Sq < j < U - 1. (13) 

1=0 ^• 
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Here, we have used the fact that qn-i,m~aj-ia„ = for all Sq < j < n — 1, 
kj + 1 <l < lo- 

Now, recall that si,0 < I < Iq is the minimal integer between 1 and n — 1 
such that m — Ian — clsi > 0, then 

n—l n—1 

Cijqn-l,m-aj-la„ = ''^^ Cijqn-l,m-aj-la„, < I < Iq. (14) 
j=so j=si 

According to (fT3l) . (fT4|) and the fact s = sq, we find that 

n—1 n—1 Iq ^ 

^ ^ CijQn,m,—aj ^ ^ Q,j ^ ^ ■^(Cl|}i'^2 ~l~ • • • ~l~ Cd^n-^d) Qn—l,m—aj—la„ 
j=s j=so 1=0 

lo ^ n-1 

= ^ ^(Cl,„X2 + . . . + Cd,nXdyC^ Cijqn-l,m-aj-laJ- 
1=0 ' j=si 



Consequently, 
dq„ 



dxi 



□ 



Next, we will give some simple examples of D-invariant subspaces as in 
Proposition [5l 

Example 1. Let Cij = 0,l<i<d,l<j<n, then 

Q = spanjrj — (ci,oXi + . . . + ^,02:^)"" : < m < ai}, 
m 

which is spanned by homogenous polynomials. We call this case the trivial 
one. □ 

Example 2. Let n = 1 and ci^ = 1, q^o = with 2 < i < d, then 

Q = spaupll, xi, ]-xl, —x1^ + ci,ixi + 02,1X2 + • • • + Cd,iXd}- 

Z CLi 

In this example, we should notice that when ai = 2, d = 2 and Ci^i = 
0, C2,i = 1, then Q = spanp{l,a;i, ^xf + X2}. This is the example that had 
been discussed by some well-known papers, for instance, [l^, p. 301], jo], p. 
811, and Illustration 6.1.9 of fH. □ 
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Example 3. Let n = 1, a = (1, 2),d = 3,ci = (1, 0), = (1, 1), Cg = (0, 1), 
then ^ ^ 

Q = spanjr{l, xi + X2, -x\ + 3:1X2 + -x\ + X2 + xs}, 
which will be considered in the next section. □ 

The next proposition not only gives us a type of interpolation point sets 
corresponding to the D-invariant subspaces as in Proposition |5], but also 
reveals their relationship. 

Proposition 6. Let a = (oq, ai, . . . , a„), Cj = (q^o, Q,!' • • • ' ^i,n) ^'^^h 1 < 
i < d, and let qn,m with < m < ai be as above, h a non-zero number in ¥. 
Then for arbitrary p e¥[x] and ^ G F"', 

0<r<m ^ ' j=0 j=0 \^^) 

iqn,m{D)pm+0{h). 

Proof. Applying Taylor Formulas, we obtain 

n n 
j=0 j=0 

c^MEi<,<d,o<,<n%7,,,£Liii^ ^^p(^) 

«:=0 ||7i||i+---+ll7d||i='= 1 " 

where = UU <7 ' 7. = (7^,o, • • • , l^,n) G N"+i, 1 < i < d. 

By means of the map r defined as in the above equality can be 
rewritten as 

n n 

p{$+{J2c,,{rhr^,...,J2^,,{rhp)) = 

j=o j=0 
^ ^ ,c7' c'fd QMi+...+\\',a\\ij)(^) 

y y {rhY \" \ II II — + o{h 

where 0(/i™'"*'^) is a polynomial in h. 



m+l ^ 
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Finally, by Lemma [H we can conclude that 



(\ n n 

u^,^„i ^ j=0 j=0 

„ „ f4^ + o(h^+'), 

r(ji,--;ld)=-m 1 d 



which leads to the proposition immediately. □ 

5. Main theorem 

In this section, we consider a particular type of ideal projectors associated 
with the above two classes of D-invariant subspaces, and constructively prove 
that C. de Boor's conjecture is true for ideal projectors of this type. 

First of all, we also introduce the notation that will be adopted in main 
theorem. 

Let ^(1), . . . , . . . , ^(^'+^) e F'^ be distinct points. For each 1 < 

k < let O^'^) C N"^ be a lower set and pC^) = {pf\ pf^) G (F"^)"^ be 
as in Section 3. Likewise, for each 1 < / < z/, let a^'^ = (cq^ a^'\ . . . , 

= (45, c2,...,c5'i(;)), 1 < z < rf, and g^'J),^, < m < af^ be as in 
Section 4. 

Theorem 7. With the notation above, let P be an ideal projector with 

ranP' = spanp{(5^(fc)D^(fe), J) : 

a G <k<fi,0<m< a'i\l <l <iy}, 

and let Ph be a Lagrange projector with 

ranP^ = spa%{(5^(fe)^^^d^^ ^^^(fe), (5^(M+o+0(O(mh) : 

a e O^'^), 1 < k < ^,0 <m < af\ 1<1 <u}, 

where h e¥\ {0} and 



i=o j=o j=o 

Then the following statements hold: 



15 



(i) There exists a positive 77 G F such that 

ranP/i = ranP, VO < \h\ < rj. 

(ii) P is the pointwise limit of Ph,0 < \h\ < r], as h tends to zero. 
Proof. Firstly, one can easily verify that 

a G 1 < k < fi,0 < m < af\l < I < u) e (¥[x]'Y 

and 

a G 1 < k < fi,0 < m < af\l < I < u) G {¥[x]'Y 
form F-bases for ranP' and ranP^ respectively, where 

fc=l k=l 

Provided that the entries of A and are arranged in the same order, namely 
for arbitrary fixed 1 < k < fi, a E D'-'^'^ the corresponding entries of A and 
A/i are in the same position, the same as for arbitrary fixed 1 < / < i^, < 
m < af\ We denote A = (Ai, . . . , A^,) and Xh = {^h,i, • • • , A^.^^), respectively. 

Let q = {qi,q2, ■ ■ ■ ,qs) be an F-basis for ranP, For convenience, we intro- 
duce two s X s matrices 

■= iKQj)l<i,j<s: ^J^Q '■= {■^h,iQj)l<i,j<s- 

and for arbitrary / G F[x], vectors 

f '■= (Aj/)i<i<s, Xj^f := {Xh,if)i<i<s- 

For arbitrary p G ¥[x], we have the following facts according to equality 
([2D and (USD. 

1. For fixed 1 < A; < /i and a G 0^^'), iD'^^k)P)i^^''^) 

can be linearly ex- 
pressed by 

d 

+ Apf ^) : /3 e 0^'^} u {0{h)} 

i=l 

since D^*^^ is lower, and moreover, the linear combination coefficient of 
each + h Yl'i=i APi'^'*) is independent of p G ¥[x]. 
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2. For fixed l</<z/andO<m< af\ can be 

linearly expressed by 

{p(|a'+0 + (f)W{rh)) : < r < m} U {0{h)}. 

Also, the linear combination coefficient of each + (j)^^\rh)) is 

independent of p G ¥[x\. 

In brief, we can conclude that there exists a nonsingular matrix T such 
that 

Kq\Kf] := T Kq|A^/] = [\^ q\X^ f] + [i5;.|e,] , (16) 

where each entry of [-E^le/i] has the same order as h. As a consequence, the 
linear systems 

{Kq) X = Kf and {Xlq) x = X^f 

are equivalent, namely they have the same set of solutions. 

(i) From f|T6|l . it follows that each entry of matrix Xj^q converges to its 
corresponding entry of matrix X^q as h tends to zero, which implies that 

linadet (a^^) = det (X^q) ■ 

Since det(A^q) 7^ 0, there exists 77 > such that 

det (Kq) 7^ 0' 0<\h\< rj. 

Notice that ffTB]) directly leads to rank (^^hq^ = rank (A^q) , 

lanPh = spanjpq, < \h\ < r], 

follows, i.e., q forms an F-basis for ranP/j. Since q is also an F-basis for ranP, 
we have 

ranP = ranP/j, < \h\ < 77. 

(ii) Suppose that x and Xq be the unique solutions of nonsingular linear 

systems 

{Xlq)x = Xlf (17) 

and 

{X^q)x = X^f (18) 
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respectively, where / G ¥[x] and < \h\ < rj. It is easy to see that 

Phf = qx and Pf = qxq. 

Notice that, as /i — ?■ 0, P is the pointwise hmit of Ph if and only if Pf is the 
coefficientwise limit of Phf for all / G F[a:;]. Therefore, it is sufficient to show 
that for every / G F[a^], the solution vector of system (fT7|l converges to the 
one of system f|T8|) as h tends to zero, namely 

lim X = xq. 

By (I16p . the linear system 



Xlq) X = XI f (19) 
can be rewritten as 

{X^q + Eh)x= {X^f + e^) . 

Since system ( IT^ is equivalent to system (fT7|l . x is also the unique solution of 
it. Consequently, using the perturbation analysis of the sensitivity of linear 
systems (see for example j25|, p. 80]), we have 



la; — ajoll < 



;A^g) 



\eh- EhXo\\+0{h^). 



Since each component of vector — E^Xq has the same order as h, it follows 
that lim 112 — XqW = 0, or, equivalently, lim 5 = Xq, which completes the 

proof of the theorem. □ 

Finally, we will give a complete example to illustrate the conclusions of 
Theorem [71 

Example 4. Let ^^^^ = (1, 1, 1), ^^2) = (q, q, 0). Let 

q = (1, X3, X2, Xi, Xg, X3X2, X3X1) 



and 



Of 

X = (5t(i),5t(i)7r ,<^f(l)7^ ''^£(1)^^ ,^f(2),6f{2){- h 7^ 

^ ^ ^ Oxi ^ 0x2 0x3 ^ ^ OXi 0X2 

S (^^\ I ^ I ^ I ^ )) 

^ 2 Oxl OxiOx2 2 0x2 0x3 
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be the F-basis for ranP and ranP', respectively. 

From Proposition [3] and Example 3, we know that this example is the 
case of Theorem [3, Therefore, we set 



1,1,1)5 C(i+/i,i,i) 



7 S{l,l+h,l), (5(1, 1,14-/1)5 ^(0,0,0)7 S{h,h'2+h,h^), '5{2fc,4fc2+2h,4h2))- 



Recalhng the proof of Theorem [71 we can obtain 
/II 1 1 1 










1 



1 





Then 



hence 



det(A^ 



h 
1 




1 



1 + h 
1 






2 + h 


7h' 



h{7h 



1 \ 

1 

1 


h) 

3) 3h / 



{h-l){2h + l){2h 



det(A^q) 7^ 0, 0<\h\< 



1 



hf 



Consequently, P is the pointwise hmit of Lagrange projector Ph, < \h\ < ^, 
as h tends to zero, with the property that ranP^ = spanj^A^. 
More precisely, we select a test function 

f{Xi, X2, X3) = 1 + (1 - + (1 - + (1 - 

to describe the perturbation procedure for the projector in this example. 
When h = 1/10,1/100,1/1000,..., we have 



^-4/ =4 



100 



34949 



439 



37867 



14520 7260 
2600449499 

1274614950^^ ~ 
722401 



X2 



9680 
46747801 

" 2549229900 
9690171 

X3X1 



Xi 



2303 
2904^ 



-X, 



233 
1452 



X2 



483294631 
121391900 



Xi 



7767 

1936 
24977753 

25492299" 



X3X1, 



1000 



25492299 2427838 
, 251000494994999 496749753001 333833081249167 

= 4 - X3 - ^^^,^^^,„^^^^^^ X2 ^„ ^1 



125249623999500 
249997752503 2 

250499247999^^ 



250499247999000 83499749333000 
747249001 667833161997 



250499247999 



XzX2 



166999498666 



Pf =4 - 2x3 - 4xi -xl + 4x3X1. 
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